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Abstract 

We study entanglement entropies of simply connected surfaces in field theories dual to 
Lovelock gravities. We consider Gauss-Bonnet and cubic Lovelock gravities in detail. In 
the conformal case the logarithmic terms in the entanglement entropy are governed by the 
conformal anomalies of the CFT; we verify that the holographic calculations are consistent 
with this property. We also compute the holographic entanglement entropy of a slab in 
the Gauss-Bonnet examples dual to relativistic and non-relativistic CFTs and discuss its 
properties. Finally, we discuss features of the entanglement entropy in the backgrounds 
dual to renormalization group flows between fixed points and comment on the implications 
for a possible c-theorem in four spacetime dimensions. 
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1. Introduction and summary 

Entanglement entropy is an interesting non-local observable which carries important infor- 
mation about field theory. Refs. [1-2] proposed a way of computing entanglement entropy 
in the strongly coupled conformal field theories dual to gravitational theories whose gravity 
sector is described by the Einstein-Hilbert lagrangian with the negative cosmological con- 
stant. The set of CFTs that admit duals of this type is strongly restricted. In particular 
in four spacetime dimensions, all such CFTs necessarily have their a and c central charges 
equal to each other. 

It is known that some interesting phenomena in CFTs and, more generally, in quantum 
field theories, are associated with the regime where a 7^ c. For example, there are unitarity 
constraints in CFTs [3] which restrict the ratio of a/c to lie within certain bounds. Another 
interesting and important question is whether one can formulate and prove the analog 
of Zamolodchikov's c-theorem [4] in three and more spacetime dimensions. It has been 
suggested that the value of the c-function in four dimensions is equal to a at fixed points 
[5]; recent work in field theory includes [6-8]. The conjecture is known by the name "a- 
theorem" . 

Holographic theories with higher derivative terms provide a natural arena for investi- 
gating these phenomena. Recently Myers and Sinha [9-10] have shown that one can formu- 
late an a-theorem in certain higher derivative theories of gravity, but it is not completely 
clear what the field theoretic counterpart of the corresponding a-function is. Their results 
have been generalized to Lovelock theories in [11]. (Work which uses higher derivative 
gravitational theories to study unitarity constraints in CFTs includes [12-25].) Interest- 
ingly, entanglement entropy provides an independent way of computing the a and c central 
charges in the CFTs. In particular, it has been noticed in [26] that depending on the shape 
of the surface which defines entanglement entropy, the logarithmic terms contain a linear 
combination of the a and c central charges. This provides an additional motivation to in- 
vestigate the holographic entanglement entropy (EE) in the theories with higher derivative 
gravitational terms. In this paper we consider Lovelock gravities, paying special attention 
to the Gauss-Bonnet and cubic Lovelock cases. We make use of the prescription of [27] 
and generalize it to the Lovelock case to compute the logarithmic terms in the holographic 
entanglement entropy for a few simple geometries such as a ball, a cylinder and a slab. 

The rest of the paper is organized as follows. In the next Section we review the results 
of [26] which imply that EE of a ball contains a logarithmic term proportional to the 
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a central charge, while the EE of a cylinder in the similar manner encodes the c central 
charge. In Section 3 we give a brief review of Lovelock theories of gravity. Section 4 contains 
the description of holographic entropy proposal of [1-2] together with the generalization 
to the Gauss-Bonnet case [27]. There we holographically compute the logarithmic terms 
in the EE of a ball and a cylinder in CFTs dual to Gauss-Bonnet gravity in AdS 5 and 
confirm that they are proportional to the a and c central charges respectively. We also 
compute the entanglement entropy of a slab as a function of Gauss-Bonnet parameter. In 
Section 5 we make an educated guess for the holographic formula valid in all Lovelock 
theories. We show that the coefficient of the logarithmic term in the EE of a cylinder in 
six-dimensional CFT dual to cubic Lovelock gravity in AdS? is proportional to the linear 
combination of the B-type anomaly coefficients. In Section 6 we consider the solution in 
the bulk which is holographically dual to the non-relativistic field theory with Lifshitz 
symmetries. We compute the entanglement entropy of a slab and a cylinder and compare 
with earlier results. Section 7 is devoted to studying entanglement entropy in the bulk 
geometries dual to renormalization group flows between conformal fixed points, motivated 
by the search for a c-theorem in four spacetime dimensions. 

Note added: As we were working on this project, we became aware of the forthcoming 
paper [28] which partially overlaps with our results (see also [29]). 

2. Entanglement entropy and conformal field theories in four dimensions. 

The entropy of entanglement (EE) in a <i-dimensional quantum field theory on IR^-i x IR 
is defined as the von Neumman entropy of the reduced density matrix associated with a 
subspace V of the total space IR^-i where the field theory lives 



The EE can be ultraviolate (UV) divergent in the continuum limit and a cutoff e needs 
to be introduced. The leading divergent term is usually proportional to the area of the 
boundary of V, (dV) 



where the proportionality coefficient 7 depends on the regularization procedure. This 
result, known as the " area law" for EE, was first found numerically [30,31] and later 



S(V) = — trypsin pv 



(2.1) 




(2.2) 
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derived analytically [32-37] . Note however that the area law is violated in the presence of 
a finite Fermi surface [38-41]. 

For (i-dimensional conformal field theories (CFTs) the structure of the divergent terms 
usually takes the following form [42] 

s(y) = 9jt0i + ... + 9jm + go[dv] lne + 8{v) (23) 

Here, s(V) is the finite part of the entropy and <7i[dV] are local, homogeneous of degree 
z, functions of the characteristic length scale of the boundary (dV). Eq. (2.3) is based 
both on the local nature of the ultraviolate divergences and on the fact that regions with 
common boundary share the same entropy. In general, the terms gd-2[dV], ■ • • , gi[dV] 
are not physical and depend on the regularization procedure. On the other hand, the 
coefficient of the logarithmic term is physical and universal in nature, not affected by 
cutoff redefinitions. 

Here our primary focus will be on conformal field theories in d = 4 dimensions where 
the universal coefficient of EE was recently obtained [26]. In particular, using the replica 
trick and conformal invariance of a four dimensional CFT on a curved manifold, [26] 
proposed that the coefficient of the logarithmically divergent term in the entanglement 
entropy of a smooth and connected region V is given by 

where (c, a) are the central charges of the four dimensional CFT and go c , go a are defined 
as follows 



goddv] = [ iV-(<<)(<<) - iV« + \r + A* 

Jov 6 
g 0a [dV] = [ R (dv) = [ R^aMn°){n]n]) - 2R^n?n» + R+ [k i k i - (k^) 2 } . 

JdV JdA 

(2.5) 

Notice that goafdV^] is simply the Euler character of the boundary manifold (dV) which we 
have also expressed (with the help of the Gauss-Codazzi identity) in terms of the ambient 
spacetime Riemann R^ uaT and Ricci R^ curvatures. (7oc[<9K| on the other hand, to the 
best of our knowledge does not have a clear geometric meaning. 

Let us clarify the notation used in (2.5). Suppose the two dimensional boundary (dV) 
is parameterized by a set of coordinates x a with a = 1, 2 whereas the spacetime metric g^ v 



where the CFT lives is spanned by coordinates X M with fj, = 0, 1, 2, 3. Then, in (2.5) 
denote two (i = 1, 2) vectors normal to the surface (#V) satisfying: 

nfnjg^ = 

8X^ v ■ ( 2 - 6 ) 

fc^„ represents the extrinsic curvature tensor of (dV) associated to the normal rii and is 
given by 

kU = -1 P X^ P < > ( 2 - 7 ) 

where 7^ represents the induced metric equal to 7^ = <7^„ — n l ^n l u . Note that v = 
0,1,2,3 are spacetime indices raised and lowered with the metric g^ v . Finally, k l is the 
trace of the extrinsic curvature tensor k % = k^g^. 

The coefficient \i in (2.5) cannot be determined by conformal invariance. In [26] 
it is fixed by requiring agreement between (2.4) and the holographic calculation. Using 
this result [26] concluded that the coefficient of the logarithmic term in the entanglement 
entropy of a ball B and a cylinder C in any four dimensional conformal field theory takes 
the form 

S(B) = ... + ^-\ne + s(B) 

yU , (2.8) 

where B, C denote a ball of radius R and an infinite cylinder of radius R and length / re- 
spectively 1 . This interesting result provides an additional characterization of the anomaly 
coefficients (c, a) through entanglement entropy. 

Independent evidence in support of the work of [26] was first given in [43] . The authors 
of [43] used Srednicki's regularization method and numerically computed the coefficient of 
the logarithmic term in the entanglement entropy of a ball for a free bosonic CFT in 
d = 4 dimensions. Their result was in complete agreement with [26]. Recently, [44] 
analytically computed the entanglement entropy for the region of a ball in a massless 
scalar field theory in arbitrary dimensions, further verifying Solodukhin's formula for this 
case. General results for the entanglement entropy for a spherical region were established 
later on in [45-48]. 

1 Note that for a region V with zero extrinsic curvature, e.g. a slab, the logarithmic term 
vanishes. 
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It is interesting to use eqs (2.4), (2.5) to compute the coefficients of the logarithmic 
term of the EE for spatial regions of different geometrical shape such as, ellipsoids, toroids 
e.t.c. It is important to stress here that the results of [26] are restricted to regions V of 
smooth geometrical shape. Otherwise, contributions from the non-smooth boundary are 
likely to modify the coefficients of the logarithmically divergent terms in the entropy [42] . 



3. Gauss-Bonnet gravity 

Among all theories of gravity which contain higher derivative terms of the Riemann tensor 
in the their action there exists a special class of theories usually referred to, as Lovelock 
gravity. This class of gravitational theories stands out both for its simplicity and the 
several properties it shares with Einstein-Hilbert gravity. In particular, it is the most 
general theory of gravity whose equations of motion involve only second order derivatives 
of the metric. It is ghost free when expanded around a Minkowski spacetime background, 
while recently, the Palatini and metric formulations of Lovelock gravity have been shown 
to be equivalent [49] . 

The action for Lovelock gravity in d + 1-dimensions is 

[-] 

S = K-T f d d+1 x^Y(-r (d ~ |£ \ P C P , (3.1) 

16ttG^ +1 J V ^ 1 (d-2)\ p p ' V ; 

where G^ 1 is the d + 1-dimensional Newton's constant, [|] denotes the integral part of |, 
X p is the p-th order Lovelock coefficient 2 and C p is the Euler density of a 2p-dimensional 
manifold. In d + 1 dimensions all C p terms with p > [|] are either total derivatives or 
vanish identically. 

In this work we will primarily focus on five dimensional gravitational theories and 
limit ourselves to the Gauss-Bonnet action. This is the simplest example of a Lovelock 
action, with only the 4-dimensional Euler density included 
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Note that in eq. (3.2) we introduced a cosmological constant term A = — and denoted 
the dimensionless Gauss-Bonnet parameter by A instead of A2 since the other Lovelock 



Note that A p are denoted as A p in [19]. 



5 



terms vanish in this case. In what follows we will retain this notation except for section 5, 
where we discuss generic Lovelock theories. The Gauss-Bonnet term £( 2 ) i n (3.2) is 



C(2) — RmnpqR MNP ® — 4RmnR MN + R 2 (3.3) 
Equations of motion derived from (3.2) are expressed in the following way 

-^9mnC + R M n + XL 2 'H% ) n = (3.4) 

with %mn defined as 

^mn = RmlpqRn LP ® — 2RmpRn P ~~ 2RmpnqR P ® + RRmn (3.5) 
Eq. (3.4) admits AdS solutions of the form[50-51] 

ds 2 = ^AdS^r + r _ m2 + j2 dx . dx i (3 6) 

AdS \ i= i J 

where the curvature scale of the AdS space is related to the cosmological constant via 3 

LAdS = — f= ol = (3.7) 

V" 1 + Vl - 4A 

Gauss-Bonnet gravity has been extensively studied in the context of the AdS/CFT 
correspondence. The basic aspects of the holographic dictionary established in the case 
of Einstein-Hilbert gravity remain the same, since the equations of motion retain their 
second order form. Moreover, the additional parameter A allows for a holographic CFT 
dual with unequal central charges (c, a). It thus provides an opportunity to investigate 
several new aspects of the correspondence (recall that all AdS backgrounds satisfying the 
Einstein-Hilbert equations of motion yield a = c). 

There are two ways to relate the gravitational parameters, the Gauss-Bonnet coupling 
A, Newton's five dimensional coupling constant and the cosmological constant L, to 
the CFT parameters (c, a). One is via a holographic calculation of the three point function 
of the stress energy tensor and the other through the holographic computation of the Weyl 
anomaly [52-53] . Both calculations yield the same result, which is a good consistency 



3 To be specific, Gauss-Bonnet gravity admits another AdS solution with a = 1 _J 1 _ AX but 
this solution is unstable and contains ghosts [50]. 
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check. The holographic calculation of the Weyl anomaly in Gauss-Bonnet gravity was 
performed in [54] . Here we simply quote the results 

c = 457r^4|f v / r^4A 



L 3 

a = 45yr^f 

Wo) 
N 



(3.8) 
2 + 3Vl - 4A] , 



where LacLS is given from (3.7). In our conventions the CFT central charges (c, a) are 
defined through the Weyl anomaly in the following way 

L 3 

It will be helpful for the calculations in the next section to have the ratio ^ff and the 
Gauss-Bonnet coefficient, A, expressed as functions of the central charges (c, a) 

L AdS _ J_(o n _ n \ x _ (a-5c)(a-c) 

Gf = ^ (iC - a) ' A= 4(a-3e)* (310) 

3c — a 

Finally, we should note that the correspondence between the positivity of the energy flux 
in a CFT [3] and causality of the boundary theory in Gauss-Bonnet gravity discussed in 
[12-15] , restricts the values of the Gauss-Bonnet parameter A to lie within the region 
— 3^3 < A < Similar results were obtained for generic Lovelock theories of gravity in 
arbitrary dimensions [16-21,25]. 



4. Holographic Entanglement Entropy Proposal 

In the context of holography, entanglement entropy received a lot of attention after the work 
of Ryu and Takayanagi where a concrete proposal for evaluating the entanglement entropy 
was set forth [1,2]. In particular, the authors of [1,2] conjectured that the entanglement 
entropy of a spatial 4 region V in a <i-dimensional CFT admitting a dual description in 
terms of Einstein-Hilbert gravity is given by 

4 Generalization to the covariant case is discussed in [55]. 



7 



where E is defined as the minimal area surface which asymptotes to the boundary of the 
spatial region V, (dV). For more details the reader is encouraged to consult [56]. 

This proposal has by now passed several tests. When for instance, the spatial region 
V extends to the whole of space, entanglement entropy should coincide with statistical 
entropy. Indeed, at finite temperature eq. (4.1) naturally reduces to the Bekenstein- 
Hawking entropy formula whereas for vanishing temperature, the dual gravitational de- 
scription contains no horizon and the entropy vanishes as it should. Other properties of 
the entanglement entropy like strong subadditivity or the fact that V and its complement 
V c have the same entropy, are also satisfied by the holographic EE formula [57]. Moreover, 
precise agreement between the holographic computation and the field theoretic one has 
been shown in the cases where explicit results are available (mostly for two dimensional 



It is interesting to generalize (4.1) to include higher derivative gravitational theories. 
The most natural idea is to replace eq. (4.1) with Wald's entropy formula. In fact, for 
Gauss-Bonnet gravity, such a proposal already exists in the literature [27] 5 . To be specific, 
the author of [27] suggested that the entanglement entropy of a connected region V of the 
dual CFT, can be computed in the case of Gauss-Bonnet gravity through the following 
formula 



Here the integral is evaluated on E, the three dimensional surface which at the boundary of 
the holographic space reduces to the two dimensional boundary (dV) of the region whose 
entropy we want to compute and which is determined by minimizing (4.2). a in the same 
expression, corresponds to the determinant of the induced metric on E whereas ife is the 
induced scalar curvature of E 6 . 

To summarize the main reasoning of [27] recall that to compute the entanglement 
entropy on the CFT side, one starts by evaluating the partition function on a <i-dimensional 
n-sheeted space - formed by gluing the n-copies of H</ along the boundary (dV). This 

5 For related work on holographic entanglement entropy and higher curvature corrections see 
also [61]. 

6 To make the variational problem well-defined a boundary term should in principle be added 
in (4.2). This term does not affect the solution of the embedding surface but it changes the value 
of the action evaluated on the solution and thus of the entanglement entropy. It turns out however 
that the boundary term only modifies the leading UV-divergent term in the entanglement entropy. 



CFTs) [58-60]. 




(4.2) 
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procedure produces a space 1Z n with conical singularities on the surface (dV) . To evaluate 
the partition function on lZ n holographically, it is necessary to identify the dual d + 1 
dimensional geometry, <S n . The latter should be a solution of the gravitational action with 
non-zero cosmological constant, which asymptotes to lZ n at the boundary. Finding such a 
solution is a difficult task. Instead, [27] assumed that <S n is given by a n-sheeted AdS^+i 
formed by gluing together n-copies of AdS^+i along a surface of codimension two. Then, 
the problem essentially reduced to that of evaluating the gravitational action functional 
on a space with conical singularities. A method for performing this calculation (at least 
in some cases) has been developed in [62] (see also [63-70]). With the use of the above 
method for Einstein-Hilbert gravity, [27] arrived at the holographic entanglement entropy 
formula of Ryu and Takayanagi. Considering Gauss-Bonnet gravity instead, leads to the 
modified expression (4.2). 

Recently [71] questioned some of the assumptions that were used in [27] to derive 
(4.2). Still, (4.2) remains a reasonable generalization of (4.1) to Gauss-Bonnet gravity. 
First, the proposal agrees with Wald's entropy formula for AdS-Schwartzchild black holes 
in Gauss-Bonnet gravity [72-75]. Therefore, whenever the spatial region V coincides with 
the total space, the entanglement entropy is guaranteed to be equal to the thermal entropy. 
Moreover, the strong subadditivity property of EE is satisfied [57]. Finally, the difference 
between (4.2) and (4.1) is the integral of a topological quantity i.e., the euler density in 
two dimensions, just like the difference between the Einstein-Hilbert and Gauss-Bonnet 
lagrangian is the Euler density term of four dimensions. 

In the following we will use the proposal of [27] to compute the entanglement entropy 
for a region bounded by ball, a cylinder and a slab. Comparison with (2.8) will provide 
yet another check of (4.2). 

4-1. The entanglement entropy of a ball. 

To compute the entanglement entropy of a ball of radius R, it is useful to parameterize 
the AdS space in the following form 



rU 2 - T 2 
as AdS — ^AdS 



+ - (~dt 2 + dr 2 + r 2 dQ 2 ) 



(4.3) 



The first step is to identify a three dimensional surface in the bulk of AdS which reduces to 
a sphere of radius R at the boundary. Taking into account the symmetries of the problem 
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we see that the surface in question is determined by a single function r(p). With this 
ansatz the induced metric of the surface can be written as follows 



rU 2 - T 2 

ab EE — n AdS 



1 



dp 2 + —dnj } . 



(4.4) 



Using (4.4) to compute the induced curvature Ry, and substituting into (4.2) yields 



s < B > = %r/V 2 



r-y/l + 4p(r'y 



1 + XaR 



(4.5) 



where a is given in (3.7) and R, the induced scalar curvature in units of the AdS radius, is 
2 [p + 4p 2 (r') 2 + 4pr (r' + 8p(r') 3 - 2pr") - r 2 (3 + 20p(r') 2 + 16pVr")] 



R = 



r 2 [l + 4p(r') 2 ] 



(4.6) 



Eq. (4.5) gives the equations of motion which determine r(p). To specify the coefficient 
of the logarithmic term it suffices to solve for r(p) to the next to leading order in the 
neighborhood of the boundary p = 0. We find that 



(4.7) 



The solution is identical in this order to the case A = 0. Substituting (4.7) into (4.5) yields 



S(B) = 



,(5) 



4G' N < Je* 



dp 



1 — 6Aa 



2p 2 



-R 1 



1 — 6Aa 
4p 



+ O(p ) 



(4.8) 



Using (3.10) and the definition of a from (3.7) we arrive at 

a R 2 



a 



(4.9) 



which is in complete agreement with (2.1 



4-2. The entanglement entropy of a cylinder. 

Consider a three dimensional surface in AdS which reduces to a two dimensional cylindrical 
surface of radius R and length / on the boundary of the AdS space. It is then natural to 
parameterize the AdS metric as follows 



d« 2 

as AdS 



T 2 

^AdS 



^! + i (- dt 2 + dz 2 + dr 2 + r 2 ^ 2 ) 
Ap 2 p v ' 



(4.10) 
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The symmetries of the problem lead us to consider a surface described by a single function 
r{p). The induced metric is then 



ds 2 - T 2 
as EE — ^AdS 



1 



l + 4p 



/ dr\ 
\P~Pj 



dp 2 + -dz 2 + r 2 $ 2 
P 



(4.11) 



Plugging this ansatz into (4.2) yields 



AG 



(5)- 

N 



2p 2 



1 + aXR 



(4.12) 



where now R is the induced curvature of the surface in units of the AdS radius is 

- _ 2 [2p(r' + 8p(r') 3 - 2pr") - r (3 + 20p(r') 2 + 16pVr")] 

r [1 + 4p(r') 2 ] 2 

The equations of motions in the vicinity of p = are satisfied by 7 



(4.13) 



Substituting (4.14) into (4.12) yields 



5(C) = ^%2,:/ / 



4G 



(5) 



(1 - 6aA).R 1 - 2a\ 
!(>/>> 



2p 2 



+ O(p ) 



With the help of (3.10) and (4.13) we arrive at 



5(C) = ^ H lne + 

V ; 90 47re 2 720 i? 



which agrees with (2.8). 



(4.14) 



(4.15) 



(4.16) 



4-3. The entanglement entropy of a slab. 

The slab geometry corresponds to the region of space bounded by— § < < § and in- 
finitely extended along the x 2 , x s directions. This is the simplest configuration to consider 
because of the large amount of symmetry. Here, it is convenient to write the AdS metric 

as 



ds z 



'AdS 



[r\ijdx % dx^ + dr 2 ) 



(4.17) 



7 Notice that the solution is again identical in this order to the case A = 0. 
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with r]ij the four dimensional Minkowski metric. The three dimensional induced surface 
can be parametrized by a single function x±(r) as follows 

ds 2 EE = ^jf- [(1 + x'^rf) dr 2 + dx\ + dx\] (4.18) 

The induced curvature of the surface (in units of the AdS radius) is non-vanishing and 
equal to 

£ = 2 3 + 3x' 1 (r) 2 + 2rx> 1 (r)x'>(r) ^ _ 

(l + x;(r) 2 ) 2 

where the primes indicate differentiation with respect to the radial coordinate. The La- 
grangian of the system is independent of x\{r) so there is a constant of motion 

A(r)(l -2 a X + ^f) = 1 
r 3 (l + xi(r) 2 ) 3 r * 

which allows us to solve for x' x {r) exactly. Since the theory is conformal, one can rescale 
the coordinate by denoting r = r/r*. Then eq. (4.20) can be written as 

-{l-2a\ + h(r)) 



(l + ftH)* 

where h(r) = x'^r) 2 . It is easy to see that there are three solutions for h(r) but only one 
of them is continuously connected with the solution of the A = case. In the following 
we restrict our attention to this solution 8 . It would be interesting to examine the other 
two solutions which at first glance appear to be complex valued. We leave this analysis to 
future work. 

We proceed to relate the constant of motion r* with the width y of the slab 

| = r* Jo (A) J (A) = J dr^hir). (4.22) 
Evaluating the action on the solution of (4.21) yields 

s(A) = %(iz^i! +7(A) i!\ (4.23) 



where 



7 (A)=4J (A) 2 f dr 
Jo 



[l + /i(r)]2 / _2aA(3 + 3/i(r) + rh'(r))\ _ (1 - 2«A) 
^ ^ [1 + h{r)] 2 J ^ 



(4.24) 



We avoid writing down the solution explicitly since it is not particularly illuminating. 
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Recall from (3.10) that A can be expressed as a function of the ratio ^ whereas A f 5 s ) 

3c— a 



4x90?r ' 



This allows us to express the final result as 



47T 90 e z y z 



(4.25) 



with 7 (a, c) = ^ g» 7(A). The result of numerical integration in eq. (4.24) is shown in 
Fig. 1. 



0.05 



^ 0.00 



-0.05 



-0.20 -0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 



Fig 1. 7(A) [see eq. (4-24)] as a function of Xa. 



We see that the EE of a slab has the form expected from field theory considerations. The 
leading divergent term follows the area law. Its coefficient is proportional to the central 
charge a just like the entanglement entropy of the ball and the infinite cylinder. However 
this coefficient depends on the regularization procedure and has no physical meaning. On 
the other hand, the coefficient of the second term in (4.25) is universal and physical. 
Computing 7(0, c) numerically we see that it changes sign approximately at aA ~ .08, 
or ^ ~ .62. The change of sign implies that 7(0, c) is not proportional to the coefficient 

ol c ~ a , and from Fig 1. 

2(Sc—a) ° 



c = (5c_q)3 °f the thermal entropy s = cT 3 . Recall also that aX 



note that 7(A) is not linear in (aX). This observation shows that 7(0, c) is not linear in a 
and c. It would be interesting to explore this further. Another direction to pursue, is to 
consider the Gauss-Bonnet holographic dual of a confining gauge theory. Computing the 
EE for a slab in this background and generalizing the results of [76] may help to clarify 
the meaning of 7 (a, c). 
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5. Generalization to Lovelock 

So far we used Fursaev's proposal to compute the leading divergent terms in the entangle- 
ment entropy of a ball, an infinite cylinder and a slab in a four dimensional CFT defined 
from Gauss-Bonnet gravity via gauge-gravity duality. Gauss-Bonnet gravity allows for the 
holographic description of CFTs with unequal central charges (c, a) giving an opportunity 
to check the proposal of [27] against the results of [26]. According to [26], the coefficient of 
the logarithmic term in the entropy should be proportional to a for a ball, c for an infinite 
cylinder and vanish for a slab, which is precisely what we found by using the prescription 
of [27]. This result led us to generalize the proposal of [27] to any Lovelock theory of 
gravity 

[-] 

S(V) = p(~r +1 (P + 1) (d (j 2 !", 2)! Vi l^C (p) . (5.1) 

Expression (5.1) coincides with the expression for the entropy of black holes in Lovelock 
gravity as established in [72-73], which in turn agrees with Wald's entropy formula [74-75]. 
As a result (5.1) satisfies by construction several properties of the entanglement entropy. 

An interesting check of the proposal (5.1) is the computation of the entanglement 
entropy of the cylinder in the case of cubic Lovelock gravity dual to the six-dimensional 
CFT. In this case, we expect the log term in the EE to be proportional to the coefficients of 
the B-type anomaly (just as in the four-dimensional case, the coefficient was proportional 
to c) 9 . In fact, there is a non-trivial check of this statement. In Lovelock gravity, the three 
coefficients of the B-type anomaly, 6j, i = 1, . . . , 3 are not linearly independent. Using the 
results of [19] one can write them as 

168 + 100A 2 a 



^AdS 


/204 




\ a 


^AdS 


(-- 




\ a 


^AdS 


(_£_ 


4^5? 


\ a 



b 2 = c^^-[ 66 + 41A 2 a (5.2) 

b 3 = c-^- (-- + 6-3X 2 a 



where we used 



\ 3 a 2 = -- + 1 - \ 2 a (5.3) 
a 



9 For the EE of a ball, the log term is shown to be proportional to the coefficient of the A-type 
anomaly in [28]. 
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to eliminate terms linear in A3. The overall coefficient c is related to the definition of the 
invariants Ii, i = 1, . . . , 3 and is therefore not important. One can fix its numerical value to 
be c = —1/2304. We can now use (5.1) to compute the EE of the cylinder. In particular, 
(4.7) becomes 

. 3p 45 - 98A 2 a - 159A 2 a 2 p 2 

r{p) = R -8R + 521(-1 + 2X ia + 3 W) W + • • • (5 ' 4) 

while the result for the EE 

94, 49 



135 L\ dS 2n 2 l 

°ios — 



4 x 512 ag (7) R 



1 - T^A 2 a - t^-A 3 q; 2 
45 15 



lne (5.5) 



Using (5.3) we can eliminate the term proportional to A 3 from the square brackets. The 
resulting expression has three terms, but turns out to be proportional to the linear com- 
bination of the two central charges (5.2), 

, 94, 49, o 19 , 1 , 

1 A 2 a A3Q; = 61 69 (5.6) 

45 15 360 10 V ; 

which is a highly nontrivial check of (5.1). Note that the third central charge, 63, is a linear 
combination of b\ and 6 2 (since Lovelock gravity was shown to satisfy the supersymmetric 
constraint [77]). Finally, the results of this section can be helpful in extending the proposal 
of [26] to CFTs of higher dimensionality. We hope to investigate this issue further in the 
future. 



6. Entanglement entropy in Lifshitz backgrounds. 

Lifshitz spacetimes are geometries of the form 

d-i 



L 

ds 2 = —2 



dr2 ~ Z2^2 df2 + J2 dX ^ 



1=1 



(6.1) 



which preserve the scaling symmetry 



r ->■ k 1 r, t-> K w t, x l ->■ kx\ . 



(6.2) 



w is usually referred to as the dynamical critical exponent. It is clear that any metric of 
the form (6.1) with w 7^ 1 breaks Lorentz invariance. 

Lifshitz spacetimes received considerable attention in the literature (together with 
Schrodinger geometries) because they provide a natural playground for the study of 
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strongly coupled systems near non-relativistic critical points [78]. Recently, Lifshitz ge- 
ometries have been embedded in string theory [79-82]. These solutions typically require 
non-trivial profile for fields other than the metric such as e.g. the dilaton. 

It is an interesting fact, that Lovelock gravity admits Lifshitz solutions of arbitrary 
critical exponent w for special values of the Lovelock parameters [83-84]. This is to be 
contrasted with the Einstein-Hilbert case, where additional matter fields are required for 
the theory to support such solutions 10 . Restricting our attention to Gauss-Bonnet gravity, 
we find that (6.1) solves the equations of motion (3.4) as long as A = \ and L 2 = 
Note that when A = \ symmetries preserved by the gravitational action functional are 
enhanced to the SO(4, 2) group. The Lagrangian of (3.2) then coincides with the Chern- 
Simons Lagrangian for the AdS group and admits a natural supersymmetrization [85]. It 
is perhaps worth mentioning that the AdS solution of Gauss-Bonnet gravity for A = \ 
is difficult to interpret in the context of the AdS/CFT correspondence since the dual 
conformal field theory has vanishing central charge c. 

Here we will focus on computing the entanglement entropy holographically in d = 4 
dimensional field theories with Lifshitz scaling. When the time coordinate is the only 
coordinate with a different scaling, the computation of entanglement entropy is identical 
to the one in AdS spacetime. Therefore the results of the previous sections go through 
unmodified as long as we make the substitution A — >■ \ and LAdS — > I n the following 
we will consider the more interesting case where rotational invariance is broken and instead 
of the time coordinate, one of the spatial coordinates e.g. x 1 , scales like x 1 — > k w x x . In 
particular, we will compute the entanglement entropy of an infinitely extended slab and a 
cylinder. Generalization to higher dimensions is straightforward. 

10 Asymptotically Lifshitz black hole solutions however do not exist in pure Lovelock gravity 
[83]. 

11 For completeness, we remark that Gauss-Bonnet gravity (for particular values of the param- 
eter A) admits Lifshitz solutions of the most generic form where several boundary coordinates 
scale like x 1 — > k w x\ However, the anisotropic scaling is a relative notion and thus, in d = 4 
there are effectively only two physically distinct cases. The first case, with a single coordinate of 
anisotropic scaling, is discussed above. The second case, of a Lifshitz spacetime with two coor- 
dinates scaling anisotropically, also satisfies the equations of motion of Gauss-Bonnet gravity as 
long as A = 1+ 1 W 2 +™ 2 and L 2 = 1+w + w2 L 2 . We will not address the latter case in the following 
which can be studied in a similar manner. 
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6.1. Entanglement entropy of an infinite belt. 

Consider double Wick rotation of (6.1), 



I?_ 

~^2 



1 / d_1 

dr2 + Z2^2 dx i+\- dt2 + T, dX 



i=2 



(6.3) 



Now there are two distinct orientations for the belt geometry, depending on whether its 
width (smallest size) extends along the direction with w-scaling or not. Here we concen- 
trate on the case where the slab is infinitely extended along the direction with anisotropic 
scaling. 

So let us consider a stripe, infinitely extended along directions x\,xs and width y 
along X2- It is convenient to parametrize the bulk surface by a single function X2(r). The 
induced metric and curvature of the surface are 



ds 2 EE 



2r 2 



R = -2 



(l + 4(rf)dr> + dxl + -^t 
(l + w + w 2 ) + (l + w + w 2 ) x' 2 (r) 2 + (1 + w)rx' 2 (r)x 2 (r) 

(i + 4M 2 ) 2 



(6.4) 



Substituting (6.4) into (4.2) results in a Lagrangian independent of X2(r) which leads to 
the following equation of motion 



x 



2 (r) (1 - w + x' 2 (r) 2 ) r 2 +* 



\2\§ 



(l + x' 2 (r) 2 ) 

Note that r* is related to the width of the slab through 



(6.5) 



- = rJ G I = J dT^/hJj) 



(6.6) 



where the integral I is expressed in terms of the dimensionless variable r = . and the 



function h(r) = (^) . 

Expressed in terms of h(r) eq. (6.5) has three solutions. Only one of them is real 
and continuously connected to the w = 1 case. Again we will restrict our attention to this 
case. Evaluating (4.2) on the particular solution of (6.5) yields the entanglement entropy 



Suf.(S) 



AG 



(5) 
N 



— W- 



h x Zc 



h x l 3 

^ yW + l 



(6.7) 
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ij. Jo is defined in (6.6) 



f 1 -w 2 - w + (1 - w - w 2 )h(r) + h(r) 2 + ^rh(r) 

/ dr 3 

Jo r 3 (l + /i(r)) 5 




(6.8) 



Let us discuss eq. (6.7). First of all, we find that the entanglement entropy is pro- 
portional to the boundary area h x l 3 of the belt. This is in accordance with expectations 
from field theory considerations [86]. We see that the logarithmically divergent term is 
absent just like in the relativistic case. The power of the leading divergent term depends on 
the critical exponent w and is fixed by dimensional analysis. The coefficient of the second 
term in (6.7) which is physical and assumed to measure the total degrees of freedom of the 
system, scales in the same way the cutoff scales. It is interesting to note that the sign of 
the leading divergent term is negative for w > in contrast to the relativistic case. This 
is not in contradiction with the results of the previous sections for Gauss-Bonnet gravity, 
because for w = 1, i.e., AdS space, the central charge c of the dual theory vanishes and a 
becomes negative (recall that Lifshitz solutions exist for A = |). For this reason the case 
w = 1 appears to be unphysical. Nevertheless, the leading divergent term in the entropy 
does not have a physical meaning and its overall sign is immaterial. We expect that a 
proper treatment of the boundary term would give a positive result but do not pursue 
this issue further here. Similar results for the entanglement entropy of a slab in a Lifshitz 
background were obtained in [87]. 

6.2. Entanglement entropy of an infinitely long cylinder. 

Another interesting case to consider is that of an infinite cylinder extended along a direction 
with anisotropic scaling. To make contact with sections 4.2 and 4.3 we express the Lifshitz 
metric in the following form 




(6.9) 



The induced metric and scalar curvature (in units -k= = 1) of the surface are 



-2(1 + w + w 2 )r(p) - 8p(2 + 2w + w 2 )r(p)r'(p) 2 



(6.10) 



r(p)(l + 4pr'(p) 2 ) 2 
Apr'{p){w + Apr'{p) 2 {l + w)) - 8p 2 r"(p)(l + 2(1 + w)r{p)r'{p)) 

r(p) (l+4pr'(p) 2 ) 2 
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where primes denote differentiation with respect to p. Combining (4.2) and (6.10) deter- 
mines the equation of motion for r(p) which takes a rather complicated form but remains 
second order in derivatives. Studying the equation of motion in the vicinity of the boundary 
p = we were able to eventually determine the exact solution. We find that 



r(p)=RJl- 



R 2 



(6.11) 



Note that (6.11) is 
for A 



independent of the critical exponent w 12 . Evaluating the action (4.2) 



j on the solution (6.11) leads to 



2~iL s 



AG 



(5) 
N 



■2nl 



I, 



, wR 
dp—j+r 
2p— 



(6.12) 



where l z regularizes the length of the cylinder. Performing the integral in (6.12) determines 
the entanglement entropy to be 



o-f r3 

S L i f .(C) = ^2irl, 



AG 



(5) 
N 



wR 

w + l 
P 2 



R 2 ) 



> w>l —2-KL, 



AG 



(5) 
N 



wR 



(6.13) 

The computation of the entanglement entropy makes sense only for w > l 13 . Recall that 
Lifshitz geometries are unstable for w < 1 (see for example [88]). We see a manifestation 
of this fact here through the computation of entanglement entropy. As expected, the 
leading divergence is proportional to the area 2nl z R whereas the scaling of the cutoff is 
fixed by dimensional analysis. Note that the logarithmic divergence, characteristic of the 
entanglement entropy on the cylinder in relativistic theories, is absent here. It would be 
interesting to examine this point further from the field theory point of view. 



12 It is interesting to expand the solution close to the boundary p = as r(p) = R — ^ — • • •. 
Note that the solution is not equal to next leading order to the embedding function for the case of 
a cylinder in Gauss-Bonnet gravity (section 4.2 eq. (4.14)). This shows once more that the order 
of limits A — >■ \ and w — > 1 cannot be interchanged. 

13 We do not consider the case w = 1 since its physical meaning is ambiguous as explained 
earlier. 
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7. Entanglement entropy and domain wall geometries 

In an attempt to investigate aspects of the entanglement entropy along renormalization 
group trajectories, we will consider here domain wall geometries. These are asymptotically 
AdS spaces with a metric of the form 

ds 2 = dr 2 + e 2A ^r ]ij dx i dx j , (7.1) 

where rjij is the Minkowski metric of the dual quantum field theory spacetime. We assume 
here that domain wall geometries are solutions of Gauss-Bonnet gravity with matter fields 
(e.g. scalars) just as they are solutions of Einstein-Hilbert gravity with matter fields. In 
the spirit of the AdS/CFT correspondence they correspond to renormalization group flow 
trajectories for the dual quantum field theories. The main difference here is that Gauss- 
Bonnet gravity contains an additional dimensionless constant A which at the fixed points 
of the flow, e.g. at r = oo, is expressed in terms of the ratio of the central charges a, c of 
the dual CFT. 

The central objective of this section is to explore how the universal and dimension- 
less coefficient of the logarithmic term in the entanglement entropy changes along the 
renormalization group flow. The main assumption we will rely on is that the holographic 
computation of entanglement entropy remains the same despite the non-trivial profile of 
the matter fields. To make the comparison with the results of sections 4.1 and 4.2 straight- 
forward, we express the metric (7.1) in a different coordinate system by performing the 
coordinate transformation p = L 2 e l . This way the IR region, r = — oo, is mapped to 
p = oo while the UV region, r = oo, corresponds to p = 0. Eq. (7.1) is expressed in the 
following form 

'dp 2 e 2U ^ 



ds 2 = L z 



H rjij dx l dx J 



_4p 

with U(p) related to A(r) in (7.1) through U(r) = A(r) — £>, or equivalently e 2A ( p > 



(7.2) 



2e 



2U(p) 



Since the metric (7.2) is asymptotically AdS, e 2U ^ admits an expansion near 



p 

,14 



the boundary p = of the form 1 

e 2U( P ) =1+(3l() + p 2 + lnp + fa ln 2 p j + . . . _ (7 3) 

The dimensionfull coefficients f3i, fa, ■ • ■ are determined by solving the field equations order 
by order in p. 
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In special cases, e.g. [89], the expansion may include half integral powers of p [90-92]. 
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7.1. The entanglement entropy of a ball along the RG flow. 



To compute the entanglement entropy of a ball of radius R we will take steps similar to 
the conformal case. Writing the domain wall metric as 



ds' 



do 2 e 2U ^ , 9 o 9 ox 

+ (-dt 2 + dr 2 + r 2 dQ 2 ) 



(7.4) 



we see that the ansatz r(p) is still natural. The induced metric is then 



EE — ^ 



Ap 2 



1+4^) (|) : 



(7.5) 



and substituting into (4.2) we arrive at 



S(B) 



L 3 fi, 



AG 



(5) 
N 



dp 



r 2 (p)e 2U ^^l + 4e 2U (p)p(r'Y 
2P 2 



1 + aXR 



(7.6) 



where R = R{r{p), U{p)) is the scalar curvature of the induced three dimensional surface 
in units of the asymptotic AdS radius in the UV, L. The equation of motion for r(p) 
derived from (7.6) is rather complicated. In the vicinity of the boundary however it is 
solved by r(p) = R — ^ + • • • exactly like the case A = 0, U(p) = 0. Knowledge of the 
near boundary behavior of r(p) together with (7.3) allows us to determine the divergent 
terms in the entanglement entropy of a ball 



S(B) 



AG 



(5) 
N 



dp 



1 — 6aA o 1 — 6aA 
R z - 



2p 2 



Ap 



1 - 2$ X R 



1 - 2aA 
1 — 6aA 



L 3 (1 - 6aA) R 2 Q 2 L 3 Q 2 1 - 6aA ( , 



AG 



(5) 



2e 2 



AG 



(5) 
N 



+ 



lne + 



(7.7) 



Here a = h-. To leading order in the vicinity of p = it is given by (3.7). 

L 2 

The coefficient of the logarithmically divergent term is modified compared to the pure 
AdS case by the overall factor ^1 — 2(3iR 2 \Zq&\ ) an d ot given in (3.7). Since ^E§^ is 
positive, the behavior of the coefficient of the logarithmic term under rescalings of the 
radius of the ball R — > AR depends on the sign of /5i . 



21 



7.2. The entanglement entropy of a cylinder along the RG flow. 

The case of a cylindrical surface is dealt with in a similar manner. We write the domain 
wall solution as 



ds 2 = L 



dp 1 
4p 



,2U(p) 



+ (-dt 2 + dz 2 + dr 2 + r 2 d<t> 2 ) 



(7- 



and select an ansatz of the form r(p) which leads to the following induced metric 



l + 4^)(|) : 



,2U(p) 



dp 2 + (dz 2 + r 2 (p)d<p 2 ) } . (7.9) 



We subsequently substitute into (4.2) and derive an equation of motion for r(p) which we 
solve near the boundary p = 0. The solution r(p) = R — ^ + (D(p 2 ), together with (7.3), 
help us obtain the divergent terms for the entanglement entropy of a cylinder 



s(c) = I dp _ (1 _ + 

AG% ] Je* L V ^pR 



L\ dS (l-6a\) 2nlR L\ dS 2nl 1 - 2a\ , 2 , 
T^Tsi ^ + 77^1 IT § (l-8/?ii2 )lne + 

N 



(7.10) 



Notice that the coefficient of the logarithmic term differs by a factor equal to 
(l — 8/3iR 2 ) compared to the conformal case. Depending on the sign of /3i, the over- 
all coefficient behaves in exactly the same way for both the ball and the cylinder. It will 
be interesting to understand the implications of this statement. 



7.3. [3i and the weak energy condition. 

Here we will examine the weak (or null) energy condition and investigate whether it is 
possible to constrain the sign of f3\ without further specifying the geometry (7.2). The null- 
energy condition implies that the matter stress energy tensor Tmn satisfies TmnC, M ' ( N > 
for any arbitrary null vector Q M . Consider the most general null vector Q M 



g^C = o =► (C p ) 2 = -A P e 2U ^ ( Vij CC j ) , 



(7.11) 



where rjij represents the four dimensional Minkowski metric. It directly follows from (7.11) 
that (vij(X j ) < 0. Substituting eq. (7.11) into the null energy condition yields 



T^C > =*► - [VijCC) (Ttt + Ape 2U ^T pp ) > . 



(7.12) 
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Here 15 we used T x i x i = —T u which follows from the symmetries of the metric (7.2). As 
we saw earlier -rjijCC > which results in (T tt + 4pe 2U ^T pp ) > 0. 

Recall that the equations of motion (3.4) for d = 4 in the presence of matter reduce 

to 

-\g MN C + Rmn + AL 2 ftgiv = (SttG^) T MN • (7.13) 
Eq. (7.13) allows the determination of (T tt + Ape 2U ^T pp ) for arbitrary matter sector 

T tt + 4pe 2U ^T pp = - 12e 2U ^ (U'(p) + pU"{p)) x 

(7.14) 

x (1 - 2\a + 8XapU'(p) [1 - pU'(p)}) 

The energy condition implies that the left hand side is positive for any p, therefore for 
p = as well. The Fefferman-Graham expansion (7.3) helps us evaluate (7.14) at p = 
and combined with the null-energy condition leads to 

a — 2 

6£i >0, (7.15) 

a 

with a = (k^j . To leading order close to the boundary a is defined by (3.7). For the 
stable AdS solution of Gauss-Bonnet gravity (a — 2) > fixes the sign of f3± to be negative. 
This is exactly what happens in Einstein-Hilbert gravity with a cosmo logical constant. It 
is interesting that had we chosen the unstable AdS solution we would have a positive 

We conclude that the coefficient of the logarithmically divergent term in the entropy of 
both the ball and the infinite cylinder is monotonically increasing with dilatations (as long 
as Pi is non- vanishing) . Unfortunately, this does not mean much. To apply the reasoning 
of [93-94] the full result, not just the logarithmic term, is required; mainly because away 
from the conformal fixed point the EE depends on two, rather than one length scale, i.e., 
the radius R of the ball and the scale of the theory p. 

7.4- Entanglement entropy of the ball and the c-theorem in Einstein-Hilbert 
gravity. 

In the previous section we saw that the coefficient of the logarithmically divergent term 
in the entanglement entropy depends only on the UV data of the theory. It is therefore 
clear that this coefficient is not a good candidate for a function decreasing along the RG 
flow and being equal to a at the fixed points. This result, together with the work of Casini 

15 Note that (7.12) can also be written as (r}ijC( j ) (-1? + T?). 
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and Huerta [42] in two-dimensional field theories, lead us to consider instead the following 
quantity 

Q = "^^fp "cft- (7.16) 

Here SEE,reg contains only the finite part of the entanglement entropy of a ball of radius R. 
We would like to examine the monotonicity properties of Q(R) holographically. A priori, 
this does not seem an easy task since we cannot compute (7.16) exactly in an arbitrary 
domain wall geometry. However, it is possible that the null energy condition and general 
characteristics of the spacetime (7.2) will determine whether (7.16) behaves monotonically 
along renormalization group trajectories. 

In the following we will attack this problem in the context of Einstein-Hilbert gravity. 
This has the advantage of being technically simpler while the relevant features of the 
problem remain the same. To simplify the analysis we will parametrize the surface in the 
bulk by p(x l ) where x % for % = 1, 2, 3 are cartesian coordinates at the boundary. This choice 
yields a lagrangian independent of x % 



S(B) = f \[dx i C{p{x i ),d i p{x i )) 

{ ° N i=1 (7.17) 



£ = 5^y / S [ ^ l)]2+4e2[/(p) ^ ) 



where p{x l ) determines the profile of the embedding surface. Translational invariance 
implies that the system has a conserved "stress-energy" tensor T- 

r -=^) (a '" ) -^ (7 ' 18) 

obeying d{T- = 0. Due to spherical symmetry T„ takes the form 

T ij = A(r)rP + 5 lJ B(r) (7.19) 
where A, B can be explicitly found to be 

z 2U{p) pl^Y e 2U{p) 1 



A( r \ 

2p 2 v /p'(r) 2 + 4e 2C/ (") / 9 2p 2 | r /( p )| ^/i + 4e 2U W pr'(p) 2 

e 2U(p) J 2U(p) J 

B(r) = —Jp>( r ) 2 +4e 2U (p)pr>(p) 2 = J I + 4e 2U (p) pr'( P y 

2p A v lp z \r {p)\ v 
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(7.20) 



In the first equality we have replaced the dependence of p on the cartesian coordinates x l 
with the spherical coordinate r. In the second equality we give A(r) , B{r) in terms of r(p) 
instead of p(r). Expressing the geometry of the embedding surface with r(p) will be more 
convenient in the following. Note that r(p) satisfies the conservation equation of T 1 - 7 , in 
other words the equations of motion derived from (7.17). The conservation equation can 
be written in a simple form in terms of the functions A(r) and B(r) 

2A(r)+r^(A(r) + B(r)) = (7.21) 
dr 

Solving for r(p) in the vicinity of the UV with the help of (7.3) and (7.20) yields 

r(p) = R ~2R + ^ + ^ bg 9 + asp3 " ' • (7 ' 22) 
where o\ is a function of the parameters (Pi, fo, A3) which specify the domain wall geometry 
but also of (T3, the coefficient of a higher order term in the near boundary expansion of 
r(p) 

_ 9 - 27fl 2 ffi + AR\2ft + 9/3 2 - 6/3 3 ) - 72R 5 a 3 

Unfortunately, <j\ cannot be determined without knowledge of the exact form of the profile 
function r(p). This will turn out to be the main difficulty in determining the behavior of 
Q(R) under rescalings of R. 

Now suppose we vary the boundary conditions in the entanglement entropy compu- 
tation as follows 

p(x l + 5x l ) = e + 5e (7.24) 

where the original boundary condition is p(x l ) = e for x % G V (here V is a sphere of radius 
R). The change in the on-shell action would then be equal to 



SS nn — 



on— shell 



4(7 l 

N 



v $dip J v 



,(5) 

with dEi the volume element on V. Substituting (7.19) into (7.25) leads to 



(7.25) 



Q = -*^ff^ == ^m R " (M R ) + B(R)) \ re9 . (7.26) 

AG N 

Having obtained the explicit form of the divergences in the previous section, it is easy to 
subtract them from A(R) + B(R) to arrive at 

Q = - R ^r = ^ml( 1 - 3R2 ^ + 8R3 ^) ( 7 - 27 ) 
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Although we have obtained a simple expression for Q(R), we cannot determine whether it 
has a monotonic behavior under rescalings of the radius R. It appears that one would need 
to know the exact solution for the profile r(p). Smoothness and other generic characteristics 
of the geometry of the induced surface did not suffice to prove the monotonicity of Q(R). 

Moreover, it is not clear how the null energy condition, pertinent to the ambient 
spacetime, will affect the behavior of Q(r) which depends on the details of the embedding 
function r(p). It is likely that even if Q(R) is a c-function in the sense of Zamolodchikov, 
it is most likely a different c-function than the standard holographic one - which in the 
coordinates used here is expressed as Choi ~ (-2t/ / (p)p+i) a ' 

Finally, the two main ingredients in the proof of [93-94] in two dimensions, are Lorentz 
invariance and the strong subadditivity property of the entanglement entropy. To make 
use of the latter, it might be more appropriate to consider the geometry of an annulus. 
We leave this issue to future investigation. 
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